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^ . ABSTRACT 

(N ■ 

, Aims. The massive twin cores NGC6334I and I(N) are in different evolutionary stages and hence ideal targets to study evolutionary variations 
C ■ within the same larger-scale environment. Here, we study the warm, compact gas components. 

■ Methods. We imaged the two regions with the Australia Telescope Compact Array (ATCA) at high angular resolution in the NH3(3,3) to (6,6) 
inversion lines. 

O I Results. Compact emission is detected toward both regions in all observed inversion lines with energy levels up to 407 K above ground. This is 
particularly surprising for NGC6334I(N) since it lacks bright infrared emission and is considered a massive cold core at an early evolutionary 
stage. High optical depth and multiply-peaked line profiles complicate rotation temperature estimates, and we can only conclude that gas com- 
ponents with temperatures > 100 K are present in both regions. Toward NGC6334I, we confirm previous reports of NH3(3,3) maser emission 



o 



^ toward the outflow bow-shocks. Furthermore, we report the first detection of an NH3(6,6) maser toward the central region of NGC6334I. 
• • This maser is centered on the second millimeter (mm) peak and elongated along the outflow axis, indicating that this mm continuum core 
_ harbors the driving source of the molecular outflow. Toward the main mm peak in NGC6334I(N), we detect a double-horn line profile in the 

, NH3(6,6) transition. The current data do not allow us to differentiate whether this double-horn profile is produced by multiple gas components 
» I along the line of sight, or whether it may trace a potential underlying massive accretion disk. The data to Figures 3 to 7 are also avail- 

^ able in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr ( 130. 79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+A/.^ 
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1. Introduction 

The massive star-forming regions NGC6334I & I(N) are lo- 
cated at the north-eastern end of the NGC6334 molecular 
cloud/Hii region complex at an approximate distance of 1 .7 kpc 



in th e southern hemisphere ( Neckel] 11978c I Straw & H viand 
19891) . The whole NGC6334 complex has been subject 
to intense studies in many wavelength bands for more 



than two decades (e.g., 


McBreen et alj 1979; Rodrieuez et al. 


1982 


; Gezaril 1982; 


Louehran et alj 


1986; dePreeetal. 


1995 


: TaDiaetal.1 19961 Sandelll 20001 Carraletal. 20021 


Beuther et al.ll2005; Hunter et al.ll2006). 
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The region NGC6334I contains the well known cometary- 
shaped Ultracompact Hii (UCHii) region NGC6334F and asso- 
ciated molecular gas and dust em ission. Various maser types 
were found, from H2O masers dMoran & Rodriguez 1980; 
Forster& Caswelll [1989) and O H masers (G aiiine^ Mute li 
19871: lBrooks& Whiteoa5l200lh to Clas s n CH3OH mase r 



emission (" Norris et alJl993tlCaswelilll997l:IWalsh et al.ll 19981) 
Furthermore, mid-infrared imaging has ide ntified a poten- 
tial exciting star for the UCHii region (iDe Buizer et al 
The main molecular gas and dust peaks are lo- 



2002) 



cated at the north-wester n edge of the UCHii region (e.g., 
Kraemer & JacksonI I1995I; iBeuther et all 120051; iHunter et al. 
20061) . Multi-frequency submillimeter (submm) and millime- 
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ter (mm) dust conti nuum imaging indicates dust temperatures 
of the order 100 K (ISandel]|l2o6ol) . Single- dish molecular line 



surveys revealed a rich line forest (e.g., 
[2OOO; Thorwirth et al. 2003; Schilk e et al. 



VlcCutcheon et al 
2006h 



compara- 
ble to those observed toward hot core regions like Orion- 
KL. A molecular outflow has been observed with a veloc- 
ity range of approximatel y 150km s ' (Bachille r & Cernicharo 



1990HLeurini et al.ll2006l) . NH3(3,3) maser as wefl as shocked 



H2 emission was reported towar d the end of th e outflow 
lobes jKr aemer & Jacksonj_ 19951: Davis & Eisloe ffel 1995; 
Persi et alJ[l996l;lMegeath & Tieftrunklll999h . 



The source NGC6334I(N) is located approximately 2' to 
the north of NGC63341. It also exhibits strong (sub)mm con- 
tinuum emission but at lower temperatures (~30K). Molecular 
line surveys of NGC63341(N) revealed that most species ex- 
hibit fainter emission than that found toward NGC63341, al- 
though there ar e some species like HC^N which are stronger in 
NGC63341(N) ('Megeath & Tieftrunk"l999^,'McCutcheon^eLal 



2OOOI; Ifhorwirth et al. 2003; Sollins & Megeath 2004). No 



mid-infrared and only weak near- infrared emission has been 
detected toward NGC6334I(N) by iTapia et all (Il996l) . Recent 



deeper near-infrared imaging h as found evidence for a cluster 
of low -mass stars in this region (iPersi et al ] l2005h . ICarral et al.l 
(l2002h report the detection of two faint cm continuum 
sources, one of them assoc iated wi th C lass 11 CH3OH maser 



emission dCaswelll 1 19971 IWalsh et al.1 (l998l) . Furthermore, 
Kogan & Slvs h' ('1998') obs erved a cluster of Class i CH3OH 



masers, and ^Megeath & Ti eftrunk ( 1999 ) detected a bipolar 
outflow in north-west south-east direction in the thermal SiO 
emission. NGC63341(N) is considered to be the younger of the 
two regions. 

To better characterize this intriguing pair of massive star- 
forming regions, we started an observational campaign with 
the Australia Telescope Compact Array (ATCA) and the 
Submillimeter Array ( SMA) from ce n timete r (cm) to submm 



and mm wavelengths. iBeuther et al.l (120051) investigated the 



regions with the ATCA in the NH3(1,1) and (2,2) lines as 
well as CH3OH emission near 25 GHz. They found com- 
pact hot gas cores in all lines toward NGC63341, however, 
in NGC63341(N) only extended NH3 emission, and no ther- 
mal CH3OH at the observed frequencies, was detected. Toward 
both regions, the temperatures were too high to derive rea- 
sonable rot ational temperatures with the low-energy-level NH3 



Hunter et al. I (l2006h identified multiple mm contin- 



transitions. 

uum sources in both regions (four in NGC63341 and seven in 
NGC63341(N)). Furthermore, an additional north-east south- 
west outflow is identified toward NGC63341(N) (Hunter et 
al. in prep.), oriented appr oximately perpendicular to the one 
previously reported by Meg eath & TieftrunkI (1 1999b . 

To get a better understanding of the hot gas components 
toward the twin cores NGC63341 and 1(N), here we report a 
follow-up study of the NH3(3,3) to (6,6) inversion lines with 
flie ATCA. 

2. Observations 

We observed NGC63341 & 1(N) in November 2005 during 
two nights with the ATCA in the compact 750D configuration. 



also including antenna 6. This results, at 25 GHz, in projected 
baselines between 3.8 and 369 k/l. The phase reference centers 
were R.A. (J2000) 17''20'"53'.44, Decl. (J2000) -35°47'02".2 
for NGC63341 and R.A. (J2000) 17''20'"54'.63, Decl. (J2000) 
-35°45'08".9 for NGC63341(N). We observed the NH3(3,3), 
(4,4), (5,5), and (6,6) inversion lines with the frequencies of 
the main hyperfine components at 23.870, 24.139, 24.533, and 
25.056GHz, respectively. The velocities relative to the local 
standard of rest (vk,-) for NGC63341 and NGC63341(N) are 

~ -7.6 and 3. 3 km s"', respectively. Good uv-coverage was 

obtained through regular switching between both sources and 
the four spectral setups. On-source integration times for each 
of the NH3(3,3) and (4,4) lines, in both sources, were 220 min- 
utes. The (5,5) and (6,6) lines were observed each for 200 min 
in NGC63341 and for 190 min in NGC63341(N). The spectral 
resolution of the observations was 62 kHz, corresponding to a 
velocity resolution of ~ 0.8 km s '. The primary beam of the 
ATCA at the given frequency is ~ 130". The data were reduced 
with the MIRIAD package. Applying a robust weighting of 1 
(closer to natural than uniform weighting, thus stressing the 
shorter baselines) the synthesized beam is 1 .9" x 1 . 1 " . The rms 
per 1 km s ' channel is ~3 mJy. 

3. Results 

The four observed NH3 inversion lines cover a broad range of 
energy levels above ground with lower energy levels Ei be- 
tween 123 and 407 K (£,(NH3(3,3) = 123 K, £,(NH3(4,4) = 
200 K, £,(NH3(5,5) = 295 K, £,( NH3(6,6) ^ 4 C I7K) The 
previous NH3(1,1) and (2,2) data of lBeuflier et al.1 (l2005l) ex- 
tend the range of observed energy levels down to 23 K. Hence, 
we are able to trace various temperature components through- 
out the two star-forming regions. All 6 inversion lines up to 
NH3(6,6) were detected in both target regions (Figs. [T] & |2]i. 
While this is less surprising forNGC6334I, the northern source 
NGC6334I(N) lacks strong infrared emission and is considered 
a prototypica l massive cold core at a very early stage of star for- 
mation (e.g.. lGezarilll982h . While it is still considered to be a 
very young source, the detection of the high energy level tran- 
sitions up to NH3(6,6) clearly shows that compact warm gas 
components already exist at this early evolutionary stage. We 
discuss both regions separately first. 



3.1. NGC6334I 

3.1 .1 . Hot thermal NH3 emission 

The uv-spectra taken on the shortest baseline in NGC6334I 
(FigID are rare examples showing even the satellite hyperfine 
components, in the highly excited NH3(5,5) and (6,6) lines, 
in emission (~26 and ~31kms"' off'set from the visr). Figure 
|3] presents integrated images of all 6 NH3 inversion transi- 
tions from the (1,1) to the (6,6) transitio n (the (1,1), (2,2) an d 
cm continuum images are adapted from iBeuther et alj|2005h . 
Except for the (3,3) line, we show the integrated emission of 
the main central hyperfine line. Only for the (3,3) transition 
we integrated over the two hyperfine components at higher fre- 
quencies (negative velocities) because the main line is contam- 
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Fig. 1. UV-Spectra of the NH3 inversion lines on the shortest baseUne toward NGC6334I. The lines mark the spectral positions 
of the main and satellite hyperfine components. 




Fig. 2. UV-Spectra of the NH3 inversion lines on the shortest baseline toward NGC6334I(N). 



inated by maser emission (see iKraemer & Jacksonlll995l and 
section [3.1.2l below). 

The NH3(1,1) and (2,2) maps show two peaks associated 
with the t wo main mm c ontinu um sources, mml and mm2, re- 
ported by iHunter et al. I (l2006h . While the two peaks have ap- 
proximately the same intensities in the low-energy NH3 lines, 
the mm dust continuum flux density of mml is appro ximately 
1.5 times higher than that of mm2 (IHunter et al.l 12006,) . The 



NH3(3,3) hyperfine map shows emission toward the north- 
eastern mm source, mml, and only weaker emission toward the 
south-western mm source, mm2. Re-inspecting the NH3(1,1) 
and (2,2) data, integrated maps of their satellite hyperfine lines 
are also stronger toward mml than toward mm2, although the 
main hyperfine lines show the same intensities toward both mm 
peaks. This indicates that the NH3 optical depth toward mm2 
has to be lower than toward mml. The higher NH3 inversion 
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lines all show a bright emission peak near mml, and a much 
weaker peak near mm2. Particularly in the NH3(4,4) and (5,5) 
lines, this weaker peak is shifted to the west of mm2 toward 
the mid-infrared source detected by De Buizer et al. (2002). 
Although in the lower energy NH3 lines the two mm peaks 
appear similar, the fact that the higher energy lines are much 
stronger toward mml indicates that this is the warmer of two 
sources. 

One of the aims of this study is to get temperature es- 
timates of the warm gas. We identify the positions of the 
two main molecular line peaks with the NH3(6,6) emission 
peak (R.A. (J2000) 17:20:53.42, Dec (J2000) -35:46:57.7, off- 
set -0.3"/4.5" from the phase center) near mml, and one 
Class n CH3OH mas er feature (R.A (J2000) 17:20:53.1 8, Dec 



Table 1. Fitted parameters and rotational temperatures in 
NGC6334I. 



Caswell 1997: Walsh et al. 1998, offset 



(J2000) -35:46:59.3, 
-3.2"/2.9" from the phase center) nearmm2. The spectra of all 
four newly observed lines toward these two positions are pre- 
sented in Fig. |4] Fitting of the hyperfine structure was done in 
CLASS which is part of the GILDAS software suiteQ. Recent 
frequency measurements for NH3 and its deuterated variants 
dCoudert & Roueflll2006 ') 1 lave been taken into account. 

Toward both positions, the NH3(3,3) and (4,4) lines are op- 
tically thick so we are not able to derive reasonable fits with 
the standard fitting procedures assuming constant excitation 
temperatures. This way we get only flat-topped synthetic spec- 
tra. Since such flat-topped spectra are not observed, this indi- 
cates that different temperatures likely exist along the line of 
sight. The most extreme case is the NH3(3,3) spectrum toward 
NGC6334I-mml (Fig. |4] top-left panel): In the opticaUy thin 
case the peak intensities of the satellite hyperfine lines should 
be 3% of the peak intensity of the main line. However, in this 
case the satellite hyperfine lines of the NH3(3,3) transition have 
more or less the same intensities as the main line, indicating the 
large optical depth. Since we cannot fit these lines well, we use 
the higher excitation NH3(5,5) and (6,6) lines for the tempera- 
ture estimates. 

Reasonable fits are possible for these two lines, and we 
measure the peak intensities, the opacities of the main lines, 
and the line-widths (Table [U. The spectral peak intensities to- 
ward these two positions are similar, which is in contrast to the 
integrated images shown in Figure [3] However, the line- widths 
toward mml are more than 2kms"' larger than those toward 
mm2; this accounts for the differences in the integrated im- 
ages. Because the NH3(5,5) and (6,6) lines are of different NH3 
species (para-NH3 and ortho-NH3, respectively), their different 
statistical weights have to be taken into account. Doing this, 
we can derive rotational temperatures (Trot) toward the molec- 
ular peak positions associated with mml and mm2. The derived 
rrot(55, 66) values are 86+20 K and 67±20K for the two posi- 
tions, respectively. Although the temperatures are comparable 
within the error-bars, the data are indicative that temperatures 
associated with mml may be higher than those associated with 
mm2. 

For the very optically thick lines where the previous fits 
failed, it is also possible to estimate approximate opacities Tmain 
from the observed line ratios between the main and the satel- 





mini 


mm 2 






- j.z/z.y 


rb(NH3(5,5)) [K] 


26.4±0.3 


32.2+0.3 


Av(NH3(5,5)) [kms-'] 


5.4±0.2 


3.3+0.2 


Vpe.k(NH3(5,5)) [kms-'] 


-7.6±0.1 


-8.1+0.1 


r„,a,„(NH3(5,5))° 


0.1±1.8 


0.4+0.5 


rb(NH3(6,6)) [K] 


15.1±0.3 


14.3+0.3 


Av(NH3(6,6)) [kms-'] 


6.4±0.6 


3.5+0.6 


Vpe.k(NH3(6,6)) [kms-'] 


-7.5±0.2 


-8.0+0.1 


r„,.in(NH3(6,6))" 


0.1±2.8 


1.4+1.6 


r„,(5,5/6, 6)* [K] 


86±20 


67+20 



" Tmain IS pfone to the largest errors, since the fits in Fig. |4] do not 
accurately reproduce the satellite lines. 

* The errors for Tio, are estimated by varying Tmai„ within errors. 



lite hyperfine components. Following Ho & TownesI ( 1983 ), we 
then calculated the rotational temperatures, Trot, and NH3 col- 
umn densities. Table|2]list the values derived via this approach. 
The estimated opacities of the (3,3) and (4,4) lines are all sig- 
nificantly larger than those we derived from the (5,5) and (6,6) 
lines. Interestingly, the rotational temperature for mm2 cor- 
responds well to the rotational temperature derived from the 
higher excited line pair above. In contrast to that, the tempera- 
ture we now derived for mml is higher than that derived from 
the higher excitation lines by a factor of about 1 .6. While the 
estimated errors associated with the the (3,3)/(4,4) hne pair are 
even larger than for the (5,5)/(6,6) pair, this derivation con- 
firms the previous assessment that mml is likely the warmer 
of the two sources. The NH3 column densities derived this way 
would correspond to H2 column densities of the order a few 
times 1 0^^ cm-^ (assuming NH3 Orio n hot core abundances as 
listed in van Dishoeck & Blake 1998h . This is about 1-2 orders 
of magnitude larger than the H2 column densities one can de- 
rive at comparable spatial scales from the mm continuum data 
( Hunter et al ] l2006h . This is an additional indicator of the lim- 
itations one has to face when the lines have as high opacities 
as in the case of NGC6334I. Future radiative transfer modeling 
may explain these observed spectral Unes better; however, this 
is beyond the scope of this paper. 



Table 2. Parameters for NGC6334I derived via the satel- 
lite/main line ratios. 



mml 



mm2 



Offset ["] 

r™ain(NH3(3,3)) 
Tmain (NH3(4,4)) 

r„,(3,3/4,4) [K] 
MNHs)" [cm-2] 



-0.3/4.5 
64+5 
28+5 

140+50 
1.0 X 10" + 30% 



-3.2/2.9 
13+5 
20+5 
70+50 
6.8 X 10"* + 30% 



See http://www.iram.fr/IRAMFR/GILDAS] 



" Since we cannot derive reasonable line-widths for these lines, for 
the column density calculations we assume Av = 5.2 and Ai' = 3.2 for 
mml and mm2, respectively (see Table[Tll. 
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Fig. 3. Integrated NH3 inversion line images of the t hermal emission tow ard NGC6334I are shown in grey-scale with solid 
contours. The (1,1) and (2,2) images are adapted from ' Beuther et al" ( 2005 ). In all cases we show the integrated main line (see 
Fig. [Til, except for the (3,3) transition where we present the emission integrated over the two hyperfine satellite components 
at blue-shifted velocities (see FigdJ. The top-row images are shown on larger spatial scales than the bottom-row images, the 
inlay size for the bottom-row images is shown in the (3,3) panel. The contouring of the (1,1) and (2,2) is done in 3o- steps 
with 3cr values of 21 and 12mJybeam"', respectively. The (3,3) to (6,6) integrated images are contoured from 20 to 90% 
(step 10%) of the peak emission with p eak values of 22.5, 32 .6, 24.2 and 17.5 mJy beam respectively. The dashed contours 
show the Ultracompact Hn region from iBeuther et al. I (l2005h . contoured from 10 to 90% (step 10%) of the peak emission of 
1 197 mJy beam The markers are labeled in the top row, SMA mm continuum em ission is from Hun ter et al. (2006), H2O 
mase r emission from [p orster & Caswelll ( 19891) . CH^O H class n maser emissi on from lWalsh et alJ (Il998i) . OH maser emission 
from Brooks & Whiteoab (l2001l) and MIR sources from toe Buizer etaL ( 2002 ). The two main mm continuum sources mml and 
mm2 are labeled in the bottom-left panel. The arrows in the top row mark the direction of the CO/SiO outflow. The synthesized 
beams are shown at the bottom-right of each panel. The data to this Figure are also available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via \http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A^ 



The NH3(6,6) line-width is larger toward mml than the 
NH3(5,5) line-width. For mm2, they are approximately the 
same within the errorbars. It would be interesting to compare 
these line- widths also with those of the lower excitation lines. 
Since, we cannot accurately fit the whole spectra (see Fig. 
we tried to derive the line widths via Gaussian fits to the satel- 
lite hyperfine components. However, for the same transitions, 
Gaussian fits of various satellite hyperfine components result 
in different line-widths (differences up to 1 kms"') indicative 
of the high optical depth and the complex source structure that 
is also evident in the various hyperfine components. Therefore, 
a comparison of the line-widths of the lower-excitation lines 



is difficult. Nevertheless, toward mml the larger line- width of 
the higher excitation NH3(6,6) line compared to NH3(5,5) indi- 
cates larger internal motions of the warmer gas which is likely 
closer to the central source. 

3.1.2. NH3(3,3) and (6,6) maser emission 

Kraemer & Jacksonid 1995b reported NH3(3,3) maser emission 
toward the outflow lobes in NGC6334I. This observation is 
confirmed by our data (Fig. |5] & |6l). The most north-eastern 
(3,3) maser spot is resolved and shows a morphology ex- 
pected from bow-shocks in molecular outflows. A Gaussian 
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fit to the spectrum extracted toward the bow-shocks results 
in a peak-brightness temperature Tf, ~ 160K and a linewidth 
Av ~ 1.6 km s"'. The maser emission is identified by a combi- 
nation of high brightness temperatures, narrow line- widths and 
peculiar spatial distributions, all with respect to the known ther- 
mal gas components. This agrees with the argument set forth 
by Kraemer & Jackson (1995) that the NH3(3,3) maser emis- 
sion is not associated with any other maser types, but that it 
is likely caused by shocks when the outflow impinges on the 
ambient molecular gas. 



Tb[K] 




Fig. 6. NH3(3,3) and (6,6) maser spectra in grey and thick lines, 
respectively. The (3,3) spectrum is extracted toward the main 
north-eastern bow-shock, and the (6,6) spectra is taken toward 
the main (6,6) maser peak close to mm2. (see Fig.|5]l. The data 
to this Figure are also available in electronic fonn at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
\http://cdsweb. u-strasbg.fr/cgi-bin/qcat ?J/A H-A/ 

Here, we report NH3(6,6) maser emission toward 
NGC6334I (Figs. H |5] & |6ll. We believe that this is the first 
detection of NH3(6,6) maser emission. The NH3(6,6) maser is 
not associated with the NH3(3,3) maser emission toward the 
outflow lobes, but rather, the (6,6) maser emission is found at 
the core center and peaks coincident with the Class ii CH3OH 
maser group associated with mm2. The maser emission is re- 
solved and elongated in the north-east south-western direction 
approximately along the large-scale outflow axis. This appar- 
ent extended structure is likely a superposition of several point- 
like sub-features which cannot be resolved by our observations. 
Weak NH3(6,6) maser emission is found in the vicinity of mml 
and associated with the H2O maser cluster. The visr of the (6,6) 
maser is +5.8 km s"' approximately 13kms"' offset from the 
thermal emission at that position. The FWHM of the maser fea- 
ture is ~ 1.3km/s, far smaller than the line-width of the ther- 
mally excited line s reported in Table | T| The H2O maser veloc- 
ities reported by jForster & Caswell [l989l) range from -45.6 
to -0.8 km s ' and are hence also offset from the NH3(6,6) 
maser velocity. Although the Class 11 CH3OH maser peaks are 
at the same positio ns, their velocities (from -6 to -12kms^', 
WalshetalJll998h are also offset from the NH3(6,6) maser 



velocity. Unfortunately, from a theoretical point of view, the 
NH3(6,6) maser has not been studied yet, and we thus have 
no predictions of what to expect. As discussed in 94.11 the 
NH3(6,6) maser morphology indicates that its emission is as- 
sociated with the molecular outflow. Although we do not yet 



understand the peculiar velocity of the (6,6) maser, this indi- 
cates that it is likely associated with some outflow-shock pro- 
cesses. Further theoretical and observational work is required 
to understand its characteristics in better detail. 



3.2. NGC6334I(N) 

One of the surprising results of our previous NH3(1,1), (2,2) 
and CH3OH (around 25 GHz) study was the absence of com- 
pact molecular line emiss ion toward the main mm continuum 
peak (iBeuther et alj|2005b . The NH3 emission of the low ex- 
citation Unes peaks approximately at the mm continuum po- 
sition, but it is extremely extended without a clearly peaked 
morphology i ndicative of compact gas cores (Fig.|2]reproduces 
the data from lBeuther et alj|2005h . Therefore, one of the aims 
of this new study was to search for compact warm gas compo- 
nents associated with the mm continuum sources, the strongest 
of which is the mml source labeled in Fig. Ill The NH3(3,3) 
emission appears less extended than the (1,1) and (2,2) maps 
but it is still relatively extended and mainly south of the mm 
continuum peaks. It is spatially associated with neither of the 
two outflows but rather located between their two southern out- 
flow lobes. Perhaps there is some kind of interaction zone be- 
tween the two molecular outflows which preferentially excites 
the NH3(3,3) line. However, this is far from clear and we re- 
frain from further interpretation of the NH3(3,3) morphology 
here. Going to the higher excitation lines, the NH3(4,4) Une is 
the first transition that shows compact emission close to the 
mm continuum peaks. However, the emission peak is offset 
by about two arcseconds to the north-west from the main mm 
continuum peak, mml. This may partially be an opacity effect 
because imaging the satellite hyperfine component at negative 
velocities - although the peak is weak on a 5cr level - it peaks 
toward the main mm source mml . Only the two highest excited 
lines, NH3(5,5) and (6,6) clearly peak toward the strongest mm 
peak mml . Both maps show a secondary peak a few arcseconds 
south-east associated with mm2. We do not detect any compact 
NH3 emission toward mm4 which is associa ted with a cm con- 
tinuum source and a Class n CH3OH maser ( Walsh et aDll998 ; 
Carral et al.l[2002HHunter et al.ll2006l) . 

Only the NH3(5,5) and (6,6) lines appear reasonable for 
a temperature estimate toward the main mm continuum and 
high-excitation NH3 peak, because the NH3(1,1) and (2,2) in- 
version lines trace only extended gas, the (3,3) line shows the 
strange emission south of the mm sources, and the (4,4) line 
peak is also offset from mml. Figure [8] presents the NH3(5,5) 
and (6,6) spectra measured toward their main emission peaks 
toward mml . While the NH3(5,5) line shows a relatively simple 
single-peaked profile, the NH3(6,6) line shows a double-peaked 
profile with the two emission peaks separated by approximately 
±2kms"' from the peak-velocity derived from the NH3(5,5) 
profile. With a Icr rms of the NH3(6,6) spectrum of 1.36 K, 
the difference between the emission peak at -1.5 km s ' and 
the dip at -3.5 km s"' is ~ 3.4cr. It is not possible to derive 
a 1 component fit for this (6,6) profile which could then be 
used with the (5,5) spectrum to derive a temperature estimate 
similar to the case in NGC6334I. Therefore, we refrain from a 
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rotational temperature estimate here, however, kinetic temper- 
atures in excess of 100 K are needed to explain the detection of 
the highly excited NH3(5,5) and (6,6) lines. 




Velocity (km/s) 



Fig. 8. NH3(5,5) and (6,6) spectra in NGC6334I(N) extracted 
after imaging toward the their peak positions associated with 
the main mm continuum peak. The spectral resolution is 
1 kms"' and the spatial resolution 2.8 x 2.4". 



4. Discussion 

4.1. The outflow driving source in NGC6334I 

What is the driving source of the large-scale north- 
east south-west outflo w in N GC6334I first reported by 
Bachiller & Cernicharo (1990) and later confirmed 
Kraemer & Jackson (1995) and Hunter et al.1 (12006.) ? 



by 



The larger line- widths observed toward mml compared to 
mm2 indicate more internal motion toward mml - either turbu- 
lent motion or organized motion due to infall, outflow or rota- 
tion. Furthermore, the fact that the highest excitation lines pref- 
erentially peak on the strongest mm continuum source mml 
indicates that this is the warmest source and potentially the 
strongest power house of the region, which would favor that 
source as the driving source of the outflow. However, mm 
line and continuum observations with the Submillimeter Array 
(Hunter et al. in prep.) show that mm2 has a higher line to con- 
tinuum ratio and that mml and mm2 are equally bright in sev- 
eral CH3OH lines with energy levels above ground of a few 
100 K. Therefore, we cannot unambiguously conclude whether 
mml is the warmer source or not. Furthermore, the NH3(6,6) 
maser peak is clearly associated with the mm continuum source 
mm2, and the position angle (PA) of the elongation of the maser 
(6,6) feature is ~48° from north, comparable to th e PA of the 
CO outflow of -46° dBachiUer & Cernicharoll990l) . This same 
orientation of the (6,6) maser elongation with the large-scale 
molecular outflow indicates that the maser emission may be 
produced by shock interaction of the outflow with the ambi- 
ent core gas close to the driving sources (in contrast to the 
NH3(3,3) maser features further outside). Therefore, we sug- 
gest that this secondary mm source, mm2, may be the driving 
source of the strong high-velocity molecular outflow. 



4.2. Double-peaked NH^(6,6) emission in 
NGC6334I(N) 

The double-peaked NH3(6,6) profile toward the main mm con- 
tinuum peak in NGC6334I(N) is peculiar and requires further 
consideration. In typical infall studies, usually an optically thin 
and an optically thick line are observed, and if the optically 
thick line shows an absorption dip with a blue-red asymme- 
try at the velocity of the peak emission of the optically thin 
line, this is a firs t indicator of potential infall motion (e.g., 
iMver s et alj 11996*). Fitting the NH3(5,5) spectrum shown in 
Figure [8] the optical depth of the main hyperfine component 
is of order unity. Although NH3(6,6) is part of the statistically 
favored ortho-NH3 (its statistical weight is double that of para- 
NH3), and we cannot precisely determine the optical depth of 
this higher excitation (6,6) line from its double-peaked profile, 
it is unlikely to be very much larger than that of the (5,5) line. 
Therefore, a typical infall interpretation appears unlikely. 

Another possibility to produce such a spectrum is the pres- 
ence of multiple gas components at different velocities, perhaps 
from an embedded binary system. In this case, the NH3(6,6) 
line has to have exactly the right optical depth to sample specif- 
ically these two gas parcels. 

Another way to get such a double-peaked profile is ro- 
tation of a central circum-protostellar disk. Such a double 
horn spectral profile is expected if one observed such a disk 
in a near-edge-on configuration where the gain path of the 
red-receding and blue-approaching disk components are long 
enough to produce enough em ission at the corresponding ve- 
locities ( Beckwith & Sargentlil99 3). The fact that we do not 
see such a double-peaked profile in the NH3(5,5) line may be 
explained by a larger optical depth of this line which hence may 
not trace the inner disk structure at all. Double-peaked profiles 
of some maser emission features were interpreted in the past 
sometime s in the fr amework of r otation in a Keplerian disk 



(e.g., Cesaronil 11990; Ponomarev et al.l 119941) . However, the 



emission observed in the NH3(6,6) line toward NGC6334I(N) 
does not appear maser-like and is likely of thermal nature. 

With the current data, we cannot distinguish whether the 
double-peaked spectrum is produced by multiple, isolated gas 
parcels or by a contiguous structure like a rotating disk. 
Nevertheless, let us follow the disk proposal for a moment: as- 
suming equilibrium between the centrifugal and gravitational 
force in an edge-on disk (inclination angle / = 90°), the radius 
of such a disk would be: 



r[AU] = 



M[Mo] 



1.13 X 10-3 Av2[kms-i] 



With a separation Av between the two NH3(6,6) peaks of 
4 km s"' and an approximate e nclosed mass of 20 estimated 
from the mm continuum data dHunter et al.ll2006l) . the approx- 
imate diameter of this tentative disk-like structure is 2200 AU, 
corresponding to an angular size of 1 .3". 

To discriminate between the diff'erent possibilities, higher 
spatial resolution observations of the NH3(6,6) line will be re- 
quired to resolve and image this potential massive accretion 
disk in more detail. Furthermore, one can look for even higher 
excitation NH3 lines: In the disk scenario, one would expect 
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line profiles similar to that observed in the (6,6) line, but with a 
larger velocity separation since the warmer gas would be found 
at smaller disk radii in a centrally heated Keplerian disk. In a 
multi-component scenario, different Une opacities could result 
in the higher energy transitions sampling diff^erent gas com- 
ponents, which would produce difi'erent line profiles for the 
higher transitions. 

5. Conclusions and Summary 

High spatial resolution imaging of the highly excited NH3 in- 
version fines (3,3) to (6,6) revealed warm gas components to- 
ward both of the twin cores NGC6334I and I(N). This is par- 
ticularly surprising for NGC63341(N) because this region was 
always considered to be the prototypical massive cold core. 
While NGC6334I(N) is still in a very young evolutionary stage 
it has already formed at least one central massive protostar that 
has begun to heat up its surrounding environment. For the well- 
known hot core NGC6334I, the detection of the higher exci- 
tation lines is less surprising. However, the optical depths of 
the NH3(1,1) to (4,4) are so high in both regions that the hy- 
perfine spectra are difficult to fit. In NGC6334I, one can use 
the (5,5) and (6,6) lines to estimate rotational temperatures 
toward the mm continuum sources mml and mm2 of 86±20 
and 67+20 K, respectively. However, since these two inversion 
lines are from different NH3 species (para- and ortho-NH3), 
such rotational temperature estimates have to been taken cau- 
tiously. We also tried to estimate the opacities, and from that the 
rotational temperatures and column densities, of the NH3(3,3) 
and (4,4) fines using the ratios between the main and satel- 
fite components of their spectra. While these estimates qual- 
itatively confirm the previous estimates they also show the 
fimitations and associated errors of the NH3 analysis in re- 
gions of that high column densities and temperatures. Toward 
NGC63341(N), a rotational temperature estimate is not possi- 
ble because the NH3(6,6) spectrum is double-peaked and hence 
does not allow reasonable line fits. Therefore, we only con- 
clude that both regions harbor dense cores with kinetic tem- 
peratures > 100 K. For better temperature estimates, one has 
either to observe optically thin isotopologues Uke '^NHs, or go 
to even higher excitation inversion transitions like the NH3(7,7) 
and NH3(8,8) lines, or try to model the spectra with more so- 
phisticated radiative transfer methods. A different approach is 
to study other molecules with different abundances and en- 
ergy levels above ground. Many such molecules have a rich 
fine spectrum in the (sub)nmi wavelength windows, and a com- 
plementary study with the SubmiUimeter Array is currently in 
progress (Hunter et al. in prep.). 

We confirm the previous detection of the NH3(3,3) maser 
toward the two outflow lobes in NGC63341. The north-eastern 
NH3(3,3) maser morphology even resembles the shape of typ- 
ical outflow bow-shocks. Furthermore, we report the first de- 
tection of NH3(6,6) maser emission toward the central cores in 
NGC6334I. The main NH3(6,6) maser peak is associated with 
the mm peak 2 and elongated along the outflow axis. Therefore, 
we suggest that the strongest line and continuum source in this 
region, mml, is not the driving source of the prominent molec- 
ular outflow, but rather the weaker sub-source, nam2. 



While past observations revealed strong line emission on 
large-spatial scales toward NGC6334I(N), no previous obser- 
vations detected compact gas cores at the center of the region. 
Observations of the highly excited NH3(4,4) to (6,6) lines for 
the first time detect compact warm gas emission from the cen- 
tral main mm peak. The observed line with the highest energy 
above ground, NH3(6,6), is not single-peaked toward that main 
mm continuum core, but it shows a double-horn profile. With 
the current data we cannot differentiate whether this double- 
horn line profile is caused by multiple gas components along 
the line of sight (maybe even a binary system) or whether 
it may be the signature of a potentially underlying massive 
accretion disk. Higher spatial resolution observations of the 
NH3(6,6) or even higher excited lines that can resolve the sub- 
structure of this core are required to solve this question. 
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- NGC6334I, mml 
Tb[K]- NH3(J.3) 



-20 
Velocity (km/s) 




Velocity (km/c) 



Fig. 4. NH3(3,3) to (6,6) spectra toward the molecular peak positions associated with the mm continuum peaks mml and mm2. 
The spectral resolution is smoothed to 1 kms"^ and the spatial resolution to 2.8" x 2.2" (the same as the spatial resolution of the 
lower energy lines presented in iBeuther et alj|2005h . The additional strong peaks at positive velocities in the NH3(6,6) spectra 
are a maser discussed in ^3.1.21 The solid lines show fits to the spectra. The (3,3) and (4,4) lines are optically thick so that 
reasonable fits are not possible. The data to this Figure are also available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr ( 130.79.128.5) or via \http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A\ 
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Fig. 5. NH3(3,3) and (6,6) maser emission in NGC6334I is shown in the left and right panel, respectively. Velocity ranges of 
integrated emission are given in each panel and the markers are the same as in Fig. [3] The contour levels are from 15 to 95% 
(step 10%) of the peak emission, which is 112.5 and 126.8 mJ y beam' ^ for the (3,3) and (6,6) maser lines. The dashed contours 
show the Ultracompact Hii region from iBeuther et al.l (l2005h . contoured from 10 to 90% (step 10%) of the peak emission of 
1 197 mJy beam The right panel again shows a smaller region as shown by the inlay in the left panel. The synthesized beams 
are shown at the bottom-right of each panel. The data to this Figure are also available in electronic form at the CDS via 



anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat? J/A+A/^ 
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Fig. 7. Integrated main hyperfine component NH3 inversion line images of the the rmal emission t oward NGC6334I(N) are shown 
in grey-scale with solid contours. The (1,1) and (2,2) images are adapted from Beuther et alj (12005 ). The top-row images are 
shown on larger spatial scales than the bottom-row images, the inlay size for the bottom-row images is shown in the (3,3) panel. 
The contour levels are 58(58)348 mJy/beam and 33(33)198 mJy/beam for the NH3(1,1) and (2,2) lines, respectively. The (3,3) 
to (6,6) integrated images are contoured from 25 to 95% (step 10%) of the peak emission with peak values of 12.9, 11.6, 10.0 
and 10.0 mJy beam \ respectively. The markers are labeled in the top row, SMA mm continuum emissio n is from Hunter e t 
al. (subm.), CH3OH class 11 maser emission from Walsh et al. ( 1 99K ) and th e cm continuum emission from lCarral et al.l (l2002h . 
The mm continuum sources mml, mm2 and mm4 from lHunter et al. I (l2006h are labeled in the bottom-left panel. The large cross 
indicates the two main axis of the two identified molecular outflows. The synthesized beams are shown at the bottom-right of 
each panel. The data to this Figure are also available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
( 130. 79.128.5) or via \http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A^ 
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^ . ABSTRACT 

(N ■ 

, Aims. The massive twin cores NGC6334I and I(N) are in different evolutionary stages and hence ideal targets to study evolutionary variations 
C ■ within the same larger-scale environment. Here, we study the warm, compact gas components. 

■ Methods. We imaged the two regions with the Australia Telescope Compact Array (ATCA) at high angular resolution in the NH3(3,3) to (6,6) 
inversion lines. 

O I Results. Compact emission is detected toward both regions in all observed inversion lines with energy levels up to 407 K above ground. This is 
particularly surprising for NGC6334I(N) since it lacks bright infrared emission and is considered a massive cold core at an early evolutionary 
stage. High optical depth and multiply-peaked line profiles complicate rotation temperature estimates, and we can only conclude that gas com- 
ponents with temperatures > 100 K are present in both regions. Toward NGC6334I, we confirm previous reports of NH3(3,3) maser emission 



o 



^ toward the outflow bow-shocks. Furthermore, we report the first detection of an NH3(6,6) maser toward the central region of NGC6334I. 
• • This maser is centered on the second millimeter (mm) peak and elongated along the outflow axis, indicating that this mm continuum core 
_ harbors the driving source of the molecular outflow. Toward the main mm peak in NGC6334I(N), we detect a double-horn line profile in the 

, NH3(6,6) transition. The current data do not allow us to differentiate whether this double-horn profile is produced by multiple gas components 
» I along the line of sight, or whether it may trace a potential underlying massive accretion disk. The data to Figures 3 to 7 are also avail- 

^ able in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr ( 130. 79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat ?J/A+A/.^ 
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1. Introduction 

The massive star-forming regions NGC6334I & I(N) are lo- 
cated at the north-eastern end of the NGC6334 molecular 
cloud/Hii region complex at an approximate distance of 1 .7 kpc 



in th e southern hemisphere ( Neckel] 11978c I Straw & H viand 
19891) . The whole NGC6334 complex has been subject 
to intense studies in many wavelength bands for more 



than two decades (e.g., 


McBreen et alj 1979; Rodrieuez et al. 


1982 


; Gezaril 1982; 


Louehran et alj 


1986; dePreeetal. 


1995 


: TaDiaetal.1 19961 Sandelll 20001 Carraletal. 20021 


Beuther et al.ll2005; Hunter et al.ll2006). 
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The region NGC6334I contains the well known cometary- 
shaped Ultracompact Hii (UCHii) region NGC6334F and asso- 
ciated molecular gas and dust em ission. Various maser types 
were found, from H2O masers dMoran & Rodriguez 1980; 
Forster& Caswelll [1989) and O H masers (G aiiine^ Mute li 
19871: lBrooks& Whiteoa5l200lh to Clas s n CH3OH mase r 



emission (" Norris et alJl993tlCaswelilll997l:IWalsh et al.ll 19981) 
Furthermore, mid-infrared imaging has ide ntified a poten- 
tial exciting star for the UCHii region (iDe Buizer et al 
The main molecular gas and dust peaks are lo- 



2002) 



cated at the north-wester n edge of the UCHii region (e.g., 
Kraemer & JacksonI I1995I; iBeuther et all 120051; iHunter et al. 
20061) . Multi-frequency submillimeter (submm) and millime- 
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ter (mm) dust conti nuum imaging indicates dust temperatures 
of the order 100 K (ISandel]|l2o6ol) . Single- dish molecular line 



surveys revealed a rich line forest (e.g., 
[2OOO; Thorwirth et al. 2003; Schilk e et al. 



VlcCutcheon et al 
2006h 



compara- 
ble to those observed toward hot core regions like Orion- 
KL. A molecular outflow has been observed with a veloc- 
ity range of approximatel y 150km s ' (Bachille r & Cernicharo 



1990HLeurini et al.ll2006l) . NH3(3,3) maser as wefl as shocked 



H2 emission was reported towar d the end of th e outflow 
lobes jKr aemer & Jacksonj_ 19951: Davis & Eisloe ffel 1995; 
Persi et alJ[l996l;lMegeath & Tieftrunklll999h . 



The source NGC6334I(N) is located approximately 2' to 
the north of NGC63341. It also exhibits strong (sub)mm con- 
tinuum emission but at lower temperatures (~30K). Molecular 
line surveys of NGC63341(N) revealed that most species ex- 
hibit fainter emission than that found toward NGC63341, al- 
though there ar e some species like HC^N which are stronger in 
NGC63341(N) ('Megeath & Tieftrunk"l999^,'McCutcheon^eLal 



2OOOI; Ifhorwirth et al. 2003; Sollins & Megeath 2004). No 



mid-infrared and only weak near- infrared emission has been 
detected toward NGC6334I(N) by iTapia et all (Il996l) . Recent 



deeper near-infrared imaging h as found evidence for a cluster 
of low -mass stars in this region (iPersi et al ] l2005h . ICarral et al.l 
(l2002h report the detection of two faint cm continuum 
sources, one of them assoc iated wi th C lass 11 CH3OH maser 



emission dCaswelll 1 19971 IWalsh et al.1 (l998l) . Furthermore, 
Kogan & Slvs h' ('1998') obs erved a cluster of Class i CH3OH 



masers, and ^Megeath & Ti eftrunk ( 1999 ) detected a bipolar 
outflow in north-west south-east direction in the thermal SiO 
emission. NGC63341(N) is considered to be the younger of the 
two regions. 

To better characterize this intriguing pair of massive star- 
forming regions, we started an observational campaign with 
the Australia Telescope Compact Array (ATCA) and the 
Submillimeter Array ( SMA) from ce n timete r (cm) to submm 



and mm wavelengths. iBeuther et al.l (120051) investigated the 



regions with the ATCA in the NH3(1,1) and (2,2) lines as 
well as CH3OH emission near 25 GHz. They found com- 
pact hot gas cores in all lines toward NGC63341, however, 
in NGC63341(N) only extended NH3 emission, and no ther- 
mal CH3OH at the observed frequencies, was detected. Toward 
both regions, the temperatures were too high to derive rea- 
sonable rot ational temperatures with the low-energy-level NH3 



Hunter et al. I (l2006h identified multiple mm contin- 



transitions. 

uum sources in both regions (four in NGC63341 and seven in 
NGC63341(N)). Furthermore, an additional north-east south- 
west outflow is identified toward NGC63341(N) (Hunter et 
al. in prep.), oriented appr oximately perpendicular to the one 
previously reported by Meg eath & TieftrunkI (1 1999b . 

To get a better understanding of the hot gas components 
toward the twin cores NGC63341 and 1(N), here we report a 
follow-up study of the NH3(3,3) to (6,6) inversion lines with 
flie ATCA. 

2. Observations 

We observed NGC63341 & 1(N) in November 2005 during 
two nights with the ATCA in the compact 750D configuration. 



also including antenna 6. This results, at 25 GHz, in projected 
baselines between 3.8 and 369 k/l. The phase reference centers 
were R.A. (J2000) 17''20'"53'.44, Decl. (J2000) -35°47'02".2 
for NGC63341 and R.A. (J2000) 17''20'"54'.63, Decl. (J2000) 
-35°45'08".9 for NGC63341(N). We observed the NH3(3,3), 
(4,4), (5,5), and (6,6) inversion lines with the frequencies of 
the main hyperfine components at 23.870, 24.139, 24.533, and 
25.056GHz, respectively. The velocities relative to the local 
standard of rest (vk,-) for NGC63341 and NGC63341(N) are 

~ -7.6 and 3. 3 km s"', respectively. Good uv-coverage was 

obtained through regular switching between both sources and 
the four spectral setups. On-source integration times for each 
of the NH3(3,3) and (4,4) lines, in both sources, were 220 min- 
utes. The (5,5) and (6,6) lines were observed each for 200 min 
in NGC63341 and for 190 min in NGC63341(N). The spectral 
resolution of the observations was 62 kHz, corresponding to a 
velocity resolution of ~ 0.8 km s '. The primary beam of the 
ATCA at the given frequency is ~ 130". The data were reduced 
with the MIRIAD package. Applying a robust weighting of 1 
(closer to natural than uniform weighting, thus stressing the 
shorter baselines) the synthesized beam is 1 .9" x 1 . 1 " . The rms 
per 1 km s ' channel is ~3 mJy. 

3. Results 

The four observed NH3 inversion lines cover a broad range of 
energy levels above ground with lower energy levels Ei be- 
tween 123 and 407 K (£,(NH3(3,3) = 123 K, £,(NH3(4,4) = 
200 K, £,(NH3(5,5) = 295 K, £,( NH3(6,6) ^ 4 C I7K) The 
previous NH3(1,1) and (2,2) data of lBeuflier et al.1 (l2005l) ex- 
tend the range of observed energy levels down to 23 K. Hence, 
we are able to trace various temperature components through- 
out the two star-forming regions. All 6 inversion lines up to 
NH3(6,6) were detected in both target regions (Figs. [T] & |2]i. 
While this is less surprising forNGC6334I, the northern source 
NGC6334I(N) lacks strong infrared emission and is considered 
a prototypica l massive cold core at a very early stage of star for- 
mation (e.g.. lGezarilll982h . While it is still considered to be a 
very young source, the detection of the high energy level tran- 
sitions up to NH3(6,6) clearly shows that compact warm gas 
components already exist at this early evolutionary stage. We 
discuss both regions separately first. 



3.1. NGC6334I 

3.1 .1 . Hot thermal NH3 emission 

The uv-spectra taken on the shortest baseline in NGC6334I 
(FigID are rare examples showing even the satellite hyperfine 
components, in the highly excited NH3(5,5) and (6,6) lines, 
in emission (~26 and ~31kms"' off'set from the visr). Figure 
|3] presents integrated images of all 6 NH3 inversion transi- 
tions from the (1,1) to the (6,6) transitio n (the (1,1), (2,2) an d 
cm continuum images are adapted from iBeuther et alj|2005h . 
Except for the (3,3) line, we show the integrated emission of 
the main central hyperfine line. Only for the (3,3) transition 
we integrated over the two hyperfine components at higher fre- 
quencies (negative velocities) because the main line is contam- 
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Fig. 1. UV-Spectra of the NH3 inversion lines on the shortest baseUne toward NGC6334I. The lines mark the spectral positions 
of the main and satellite hyperfine components. 




Fig. 2. UV-Spectra of the NH3 inversion lines on the shortest baseline toward NGC6334I(N). 



inated by maser emission (see iKraemer & Jacksonlll995l and 
section [3.1.2l below). 

The NH3(1,1) and (2,2) maps show two peaks associated 
with the t wo main mm c ontinu um sources, mml and mm2, re- 
ported by iHunter et al. I (l2006h . While the two peaks have ap- 
proximately the same intensities in the low-energy NH3 lines, 
the mm dust continuum flux density of mml is appro ximately 
1.5 times higher than that of mm2 (IHunter et al.l 12006,) . The 



NH3(3,3) hyperfine map shows emission toward the north- 
eastern mm source, mml, and only weaker emission toward the 
south-western mm source, mm2. Re-inspecting the NH3(1,1) 
and (2,2) data, integrated maps of their satellite hyperfine lines 
are also stronger toward mml than toward mm2, although the 
main hyperfine lines show the same intensities toward both mm 
peaks. This indicates that the NH3 optical depth toward mm2 
has to be lower than toward mml. The higher NH3 inversion 
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lines all show a bright emission peak near mml, and a much 
weaker peak near mm2. Particularly in the NH3(4,4) and (5,5) 
lines, this weaker peak is shifted to the west of mm2 toward 
the mid-infrared source detected by De Buizer et al. (2002). 
Although in the lower energy NH3 lines the two mm peaks 
appear similar, the fact that the higher energy lines are much 
stronger toward mml indicates that this is the warmer of two 
sources. 

One of the aims of this study is to get temperature es- 
timates of the warm gas. We identify the positions of the 
two main molecular line peaks with the NH3(6,6) emission 
peak (R.A. (J2000) 17:20:53.42, Dec (J2000) -35:46:57.7, off- 
set -0.3"/4.5" from the phase center) near mml, and one 
Class n CH3OH mas er feature (R.A (J2000) 17:20:53.1 8, Dec 



Table 1. Fitted parameters and rotational temperatures in 
NGC6334I. 



Caswell 1997: Walsh et al. 1998, offset 



(J2000) -35:46:59.3, 
-3.2"/2.9" from the phase center) nearmm2. The spectra of all 
four newly observed lines toward these two positions are pre- 
sented in Fig. |4] Fitting of the hyperfine structure was done in 
CLASS which is part of the GILDAS software suiteQ. Recent 
frequency measurements for NH3 and its deuterated variants 
dCoudert & Roueflll2006 ') 1 lave been taken into account. 

Toward both positions, the NH3(3,3) and (4,4) lines are op- 
tically thick so we are not able to derive reasonable fits with 
the standard fitting procedures assuming constant excitation 
temperatures. This way we get only flat-topped synthetic spec- 
tra. Since such flat-topped spectra are not observed, this indi- 
cates that different temperatures likely exist along the line of 
sight. The most extreme case is the NH3(3,3) spectrum toward 
NGC6334I-mml (Fig. |4] top-left panel): In the opticaUy thin 
case the peak intensities of the satellite hyperfine lines should 
be 3% of the peak intensity of the main line. However, in this 
case the satellite hyperfine lines of the NH3(3,3) transition have 
more or less the same intensities as the main line, indicating the 
large optical depth. Since we cannot fit these lines well, we use 
the higher excitation NH3(5,5) and (6,6) lines for the tempera- 
ture estimates. 

Reasonable fits are possible for these two lines, and we 
measure the peak intensities, the opacities of the main lines, 
and the line-widths (Table [U. The spectral peak intensities to- 
ward these two positions are similar, which is in contrast to the 
integrated images shown in Figure [3] However, the line- widths 
toward mml are more than 2kms"' larger than those toward 
mm2; this accounts for the differences in the integrated im- 
ages. Because the NH3(5,5) and (6,6) lines are of different NH3 
species (para-NH3 and ortho-NH3, respectively), their different 
statistical weights have to be taken into account. Doing this, 
we can derive rotational temperatures (Trot) toward the molec- 
ular peak positions associated with mml and mm2. The derived 
rrot(55, 66) values are 86+20 K and 67±20K for the two posi- 
tions, respectively. Although the temperatures are comparable 
within the error-bars, the data are indicative that temperatures 
associated with mml may be higher than those associated with 
mm2. 

For the very optically thick lines where the previous fits 
failed, it is also possible to estimate approximate opacities Tmain 
from the observed line ratios between the main and the satel- 





mini 


mm 2 






- j.z/z.y 


rb(NH3(5,5)) [K] 


26.4±0.3 


32.2+0.3 


Av(NH3(5,5)) [kms-'] 


5.4±0.2 


3.3+0.2 


Vpe.k(NH3(5,5)) [kms-'] 


-7.6±0.1 


-8.1+0.1 


r„,a,„(NH3(5,5))° 


0.1±1.8 


0.4+0.5 


rb(NH3(6,6)) [K] 


15.1±0.3 


14.3+0.3 


Av(NH3(6,6)) [kms-'] 


6.4±0.6 


3.5+0.6 


Vpe.k(NH3(6,6)) [kms-'] 


-7.5±0.2 


-8.0+0.1 


r„,.in(NH3(6,6))" 


0.1±2.8 


1.4+1.6 


r„,(5,5/6, 6)* [K] 


86±20 


67+20 



" Tmain IS pfone to the largest errors, since the fits in Fig. |4] do not 
accurately reproduce the satellite lines. 

* The errors for Tio, are estimated by varying Tmai„ within errors. 



lite hyperfine components. Following Ho & TownesI ( 1983 ), we 
then calculated the rotational temperatures, Trot, and NH3 col- 
umn densities. Table|2]list the values derived via this approach. 
The estimated opacities of the (3,3) and (4,4) lines are all sig- 
nificantly larger than those we derived from the (5,5) and (6,6) 
lines. Interestingly, the rotational temperature for mm2 cor- 
responds well to the rotational temperature derived from the 
higher excited line pair above. In contrast to that, the tempera- 
ture we now derived for mml is higher than that derived from 
the higher excitation lines by a factor of about 1 .6. While the 
estimated errors associated with the the (3,3)/(4,4) hne pair are 
even larger than for the (5,5)/(6,6) pair, this derivation con- 
firms the previous assessment that mml is likely the warmer 
of the two sources. The NH3 column densities derived this way 
would correspond to H2 column densities of the order a few 
times 1 0^^ cm-^ (assuming NH3 Orio n hot core abundances as 
listed in van Dishoeck & Blake 1998h . This is about 1-2 orders 
of magnitude larger than the H2 column densities one can de- 
rive at comparable spatial scales from the mm continuum data 
( Hunter et al ] l2006h . This is an additional indicator of the lim- 
itations one has to face when the lines have as high opacities 
as in the case of NGC6334I. Future radiative transfer modeling 
may explain these observed spectral Unes better; however, this 
is beyond the scope of this paper. 



Table 2. Parameters for NGC6334I derived via the satel- 
lite/main line ratios. 



mml 



mm2 



Offset ["] 

r™ain(NH3(3,3)) 
Tmain (NH3(4,4)) 

r„,(3,3/4,4) [K] 
MNHs)" [cm-2] 



-0.3/4.5 
64+5 
28+5 

140+50 
1.0 X 10" + 30% 



-3.2/2.9 
13+5 
20+5 
70+50 
6.8 X 10"* + 30% 



See http://www.iram.fr/IRAMFR/GILDAS] 



" Since we cannot derive reasonable line-widths for these lines, for 
the column density calculations we assume Av = 5.2 and Ai' = 3.2 for 
mml and mm2, respectively (see Table[Tll. 
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Fig. 3. Integrated NH3 inversion line images of the t hermal emission tow ard NGC6334I are shown in grey-scale with solid 
contours. The (1,1) and (2,2) images are adapted from ' Beuther et al" ( 2005 ). In all cases we show the integrated main line (see 
Fig. [Til, except for the (3,3) transition where we present the emission integrated over the two hyperfine satellite components 
at blue-shifted velocities (see FigdJ. The top-row images are shown on larger spatial scales than the bottom-row images, the 
inlay size for the bottom-row images is shown in the (3,3) panel. The contouring of the (1,1) and (2,2) is done in 3o- steps 
with 3cr values of 21 and 12mJybeam"', respectively. The (3,3) to (6,6) integrated images are contoured from 20 to 90% 
(step 10%) of the peak emission with p eak values of 22.5, 32 .6, 24.2 and 17.5 mJy beam respectively. The dashed contours 
show the Ultracompact Hn region from iBeuther et al. I (l2005h . contoured from 10 to 90% (step 10%) of the peak emission of 
1 197 mJy beam The markers are labeled in the top row, SMA mm continuum em ission is from Hun ter et al. (2006), H2O 
mase r emission from [p orster & Caswelll ( 19891) . CH^O H class n maser emissi on from lWalsh et alJ (Il998i) . OH maser emission 
from Brooks & Whiteoab (l2001l) and MIR sources from toe Buizer etaL ( 2002 ). The two main mm continuum sources mml and 
mm2 are labeled in the bottom-left panel. The arrows in the top row mark the direction of the CO/SiO outflow. The synthesized 
beams are shown at the bottom-right of each panel. The data to this Figure are also available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via \http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A^ 



The NH3(6,6) line-width is larger toward mml than the 
NH3(5,5) line-width. For mm2, they are approximately the 
same within the errorbars. It would be interesting to compare 
these line- widths also with those of the lower excitation lines. 
Since, we cannot accurately fit the whole spectra (see Fig. 
we tried to derive the line widths via Gaussian fits to the satel- 
lite hyperfine components. However, for the same transitions, 
Gaussian fits of various satellite hyperfine components result 
in different line-widths (differences up to 1 kms"') indicative 
of the high optical depth and the complex source structure that 
is also evident in the various hyperfine components. Therefore, 
a comparison of the line-widths of the lower-excitation lines 



is difficult. Nevertheless, toward mml the larger line- width of 
the higher excitation NH3(6,6) line compared to NH3(5,5) indi- 
cates larger internal motions of the warmer gas which is likely 
closer to the central source. 

3.1.2. NH3(3,3) and (6,6) maser emission 

Kraemer & Jacksonid 1995b reported NH3(3,3) maser emission 
toward the outflow lobes in NGC6334I. This observation is 
confirmed by our data (Fig. |5] & |6l). The most north-eastern 
(3,3) maser spot is resolved and shows a morphology ex- 
pected from bow-shocks in molecular outflows. A Gaussian 
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fit to the spectrum extracted toward the bow-shocks results 
in a peak-brightness temperature Tf, ~ 160K and a linewidth 
Av ~ 1.6 km s"'. The maser emission is identified by a combi- 
nation of high brightness temperatures, narrow line- widths and 
peculiar spatial distributions, all with respect to the known ther- 
mal gas components. This agrees with the argument set forth 
by Kraemer & Jackson (1995) that the NH3(3,3) maser emis- 
sion is not associated with any other maser types, but that it 
is likely caused by shocks when the outflow impinges on the 
ambient molecular gas. 



Tb[K] 




Fig. 6. NH3(3,3) and (6,6) maser spectra in grey and thick lines, 
respectively. The (3,3) spectrum is extracted toward the main 
north-eastern bow-shock, and the (6,6) spectra is taken toward 
the main (6,6) maser peak close to mm2. (see Fig.|5]l. The data 
to this Figure are also available in electronic fonn at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
\http://cdsweb. u-strasbg.fr/cgi-bin/qcat ?J/A H-A/ 

Here, we report NH3(6,6) maser emission toward 
NGC6334I (Figs. H |5] & |6ll. We believe that this is the first 
detection of NH3(6,6) maser emission. The NH3(6,6) maser is 
not associated with the NH3(3,3) maser emission toward the 
outflow lobes, but rather, the (6,6) maser emission is found at 
the core center and peaks coincident with the Class ii CH3OH 
maser group associated with mm2. The maser emission is re- 
solved and elongated in the north-east south-western direction 
approximately along the large-scale outflow axis. This appar- 
ent extended structure is likely a superposition of several point- 
like sub-features which cannot be resolved by our observations. 
Weak NH3(6,6) maser emission is found in the vicinity of mml 
and associated with the H2O maser cluster. The visr of the (6,6) 
maser is +5.8 km s"' approximately 13kms"' offset from the 
thermal emission at that position. The FWHM of the maser fea- 
ture is ~ 1.3km/s, far smaller than the line-width of the ther- 
mally excited line s reported in Table | T| The H2O maser veloc- 
ities reported by jForster & Caswell [l989l) range from -45.6 
to -0.8 km s ' and are hence also offset from the NH3(6,6) 
maser velocity. Although the Class 11 CH3OH maser peaks are 
at the same positio ns, their velocities (from -6 to -12kms^', 
WalshetalJll998h are also offset from the NH3(6,6) maser 



velocity. Unfortunately, from a theoretical point of view, the 
NH3(6,6) maser has not been studied yet, and we thus have 
no predictions of what to expect. As discussed in 94.11 the 
NH3(6,6) maser morphology indicates that its emission is as- 
sociated with the molecular outflow. Although we do not yet 



understand the peculiar velocity of the (6,6) maser, this indi- 
cates that it is likely associated with some outflow-shock pro- 
cesses. Further theoretical and observational work is required 
to understand its characteristics in better detail. 



3.2. NGC6334I(N) 

One of the surprising results of our previous NH3(1,1), (2,2) 
and CH3OH (around 25 GHz) study was the absence of com- 
pact molecular line emiss ion toward the main mm continuum 
peak (iBeuther et alj|2005b . The NH3 emission of the low ex- 
citation Unes peaks approximately at the mm continuum po- 
sition, but it is extremely extended without a clearly peaked 
morphology i ndicative of compact gas cores (Fig.|2]reproduces 
the data from lBeuther et alj|2005h . Therefore, one of the aims 
of this new study was to search for compact warm gas compo- 
nents associated with the mm continuum sources, the strongest 
of which is the mml source labeled in Fig. Ill The NH3(3,3) 
emission appears less extended than the (1,1) and (2,2) maps 
but it is still relatively extended and mainly south of the mm 
continuum peaks. It is spatially associated with neither of the 
two outflows but rather located between their two southern out- 
flow lobes. Perhaps there is some kind of interaction zone be- 
tween the two molecular outflows which preferentially excites 
the NH3(3,3) line. However, this is far from clear and we re- 
frain from further interpretation of the NH3(3,3) morphology 
here. Going to the higher excitation lines, the NH3(4,4) Une is 
the first transition that shows compact emission close to the 
mm continuum peaks. However, the emission peak is offset 
by about two arcseconds to the north-west from the main mm 
continuum peak, mml. This may partially be an opacity effect 
because imaging the satellite hyperfine component at negative 
velocities - although the peak is weak on a 5cr level - it peaks 
toward the main mm source mml . Only the two highest excited 
lines, NH3(5,5) and (6,6) clearly peak toward the strongest mm 
peak mml . Both maps show a secondary peak a few arcseconds 
south-east associated with mm2. We do not detect any compact 
NH3 emission toward mm4 which is associa ted with a cm con- 
tinuum source and a Class n CH3OH maser ( Walsh et aDll998 ; 
Carral et al.l[2002HHunter et al.ll2006l) . 

Only the NH3(5,5) and (6,6) lines appear reasonable for 
a temperature estimate toward the main mm continuum and 
high-excitation NH3 peak, because the NH3(1,1) and (2,2) in- 
version lines trace only extended gas, the (3,3) line shows the 
strange emission south of the mm sources, and the (4,4) line 
peak is also offset from mml. Figure [8] presents the NH3(5,5) 
and (6,6) spectra measured toward their main emission peaks 
toward mml . While the NH3(5,5) line shows a relatively simple 
single-peaked profile, the NH3(6,6) line shows a double-peaked 
profile with the two emission peaks separated by approximately 
±2kms"' from the peak-velocity derived from the NH3(5,5) 
profile. With a Icr rms of the NH3(6,6) spectrum of 1.36 K, 
the difference between the emission peak at -1.5 km s ' and 
the dip at -3.5 km s"' is ~ 3.4cr. It is not possible to derive 
a 1 component fit for this (6,6) profile which could then be 
used with the (5,5) spectrum to derive a temperature estimate 
similar to the case in NGC6334I. Therefore, we refrain from a 
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rotational temperature estimate here, however, kinetic temper- 
atures in excess of 100 K are needed to explain the detection of 
the highly excited NH3(5,5) and (6,6) lines. 




Velocity (km/s) 



Fig. 8. NH3(5,5) and (6,6) spectra in NGC6334I(N) extracted 
after imaging toward the their peak positions associated with 
the main mm continuum peak. The spectral resolution is 
1 kms"' and the spatial resolution 2.8 x 2.4". 



4. Discussion 

4.1. The outflow driving source in NGC6334I 

What is the driving source of the large-scale north- 
east south-west outflo w in N GC6334I first reported by 
Bachiller & Cernicharo (1990) and later confirmed 
Kraemer & Jackson (1995) and Hunter et al.1 (12006.) ? 



by 



The larger line- widths observed toward mml compared to 
mm2 indicate more internal motion toward mml - either turbu- 
lent motion or organized motion due to infall, outflow or rota- 
tion. Furthermore, the fact that the highest excitation lines pref- 
erentially peak on the strongest mm continuum source mml 
indicates that this is the warmest source and potentially the 
strongest power house of the region, which would favor that 
source as the driving source of the outflow. However, mm 
line and continuum observations with the Submillimeter Array 
(Hunter et al. in prep.) show that mm2 has a higher line to con- 
tinuum ratio and that mml and mm2 are equally bright in sev- 
eral CH3OH lines with energy levels above ground of a few 
100 K. Therefore, we cannot unambiguously conclude whether 
mml is the warmer source or not. Furthermore, the NH3(6,6) 
maser peak is clearly associated with the mm continuum source 
mm2, and the position angle (PA) of the elongation of the maser 
(6,6) feature is ~48° from north, comparable to th e PA of the 
CO outflow of -46° dBachiUer & Cernicharoll990l) . This same 
orientation of the (6,6) maser elongation with the large-scale 
molecular outflow indicates that the maser emission may be 
produced by shock interaction of the outflow with the ambi- 
ent core gas close to the driving sources (in contrast to the 
NH3(3,3) maser features further outside). Therefore, we sug- 
gest that this secondary mm source, mm2, may be the driving 
source of the strong high-velocity molecular outflow. 



4.2. Double-peaked NH^(6,6) emission in 
NGC6334I(N) 

The double-peaked NH3(6,6) profile toward the main mm con- 
tinuum peak in NGC6334I(N) is peculiar and requires further 
consideration. In typical infall studies, usually an optically thin 
and an optically thick line are observed, and if the optically 
thick line shows an absorption dip with a blue-red asymme- 
try at the velocity of the peak emission of the optically thin 
line, this is a firs t indicator of potential infall motion (e.g., 
iMver s et alj 11996*). Fitting the NH3(5,5) spectrum shown in 
Figure [8] the optical depth of the main hyperfine component 
is of order unity. Although NH3(6,6) is part of the statistically 
favored ortho-NH3 (its statistical weight is double that of para- 
NH3), and we cannot precisely determine the optical depth of 
this higher excitation (6,6) line from its double-peaked profile, 
it is unlikely to be very much larger than that of the (5,5) line. 
Therefore, a typical infall interpretation appears unlikely. 

Another possibility to produce such a spectrum is the pres- 
ence of multiple gas components at different velocities, perhaps 
from an embedded binary system. In this case, the NH3(6,6) 
line has to have exactly the right optical depth to sample specif- 
ically these two gas parcels. 

Another way to get such a double-peaked profile is ro- 
tation of a central circum-protostellar disk. Such a double 
horn spectral profile is expected if one observed such a disk 
in a near-edge-on configuration where the gain path of the 
red-receding and blue-approaching disk components are long 
enough to produce enough em ission at the corresponding ve- 
locities ( Beckwith & Sargentlil99 3). The fact that we do not 
see such a double-peaked profile in the NH3(5,5) line may be 
explained by a larger optical depth of this line which hence may 
not trace the inner disk structure at all. Double-peaked profiles 
of some maser emission features were interpreted in the past 
sometime s in the fr amework of r otation in a Keplerian disk 



(e.g., Cesaronil 11990; Ponomarev et al.l 119941) . However, the 



emission observed in the NH3(6,6) line toward NGC6334I(N) 
does not appear maser-like and is likely of thermal nature. 

With the current data, we cannot distinguish whether the 
double-peaked spectrum is produced by multiple, isolated gas 
parcels or by a contiguous structure like a rotating disk. 
Nevertheless, let us follow the disk proposal for a moment: as- 
suming equilibrium between the centrifugal and gravitational 
force in an edge-on disk (inclination angle / = 90°), the radius 
of such a disk would be: 



r[AU] = 



M[Mo] 



1.13 X 10-3 Av2[kms-i] 



With a separation Av between the two NH3(6,6) peaks of 
4 km s"' and an approximate e nclosed mass of 20 estimated 
from the mm continuum data dHunter et al.ll2006l) . the approx- 
imate diameter of this tentative disk-like structure is 2200 AU, 
corresponding to an angular size of 1 .3". 

To discriminate between the diff'erent possibilities, higher 
spatial resolution observations of the NH3(6,6) line will be re- 
quired to resolve and image this potential massive accretion 
disk in more detail. Furthermore, one can look for even higher 
excitation NH3 lines: In the disk scenario, one would expect 
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line profiles similar to that observed in the (6,6) line, but with a 
larger velocity separation since the warmer gas would be found 
at smaller disk radii in a centrally heated Keplerian disk. In a 
multi-component scenario, different Une opacities could result 
in the higher energy transitions sampling diff^erent gas com- 
ponents, which would produce difi'erent line profiles for the 
higher transitions. 

5. Conclusions and Summary 

High spatial resolution imaging of the highly excited NH3 in- 
version fines (3,3) to (6,6) revealed warm gas components to- 
ward both of the twin cores NGC6334I and I(N). This is par- 
ticularly surprising for NGC63341(N) because this region was 
always considered to be the prototypical massive cold core. 
While NGC6334I(N) is still in a very young evolutionary stage 
it has already formed at least one central massive protostar that 
has begun to heat up its surrounding environment. For the well- 
known hot core NGC6334I, the detection of the higher exci- 
tation lines is less surprising. However, the optical depths of 
the NH3(1,1) to (4,4) are so high in both regions that the hy- 
perfine spectra are difficult to fit. In NGC6334I, one can use 
the (5,5) and (6,6) lines to estimate rotational temperatures 
toward the mm continuum sources mml and mm2 of 86±20 
and 67+20 K, respectively. However, since these two inversion 
lines are from different NH3 species (para- and ortho-NH3), 
such rotational temperature estimates have to been taken cau- 
tiously. We also tried to estimate the opacities, and from that the 
rotational temperatures and column densities, of the NH3(3,3) 
and (4,4) fines using the ratios between the main and satel- 
fite components of their spectra. While these estimates qual- 
itatively confirm the previous estimates they also show the 
fimitations and associated errors of the NH3 analysis in re- 
gions of that high column densities and temperatures. Toward 
NGC63341(N), a rotational temperature estimate is not possi- 
ble because the NH3(6,6) spectrum is double-peaked and hence 
does not allow reasonable line fits. Therefore, we only con- 
clude that both regions harbor dense cores with kinetic tem- 
peratures > 100 K. For better temperature estimates, one has 
either to observe optically thin isotopologues Uke '^NHs, or go 
to even higher excitation inversion transitions like the NH3(7,7) 
and NH3(8,8) lines, or try to model the spectra with more so- 
phisticated radiative transfer methods. A different approach is 
to study other molecules with different abundances and en- 
ergy levels above ground. Many such molecules have a rich 
fine spectrum in the (sub)nmi wavelength windows, and a com- 
plementary study with the SubmiUimeter Array is currently in 
progress (Hunter et al. in prep.). 

We confirm the previous detection of the NH3(3,3) maser 
toward the two outflow lobes in NGC63341. The north-eastern 
NH3(3,3) maser morphology even resembles the shape of typ- 
ical outflow bow-shocks. Furthermore, we report the first de- 
tection of NH3(6,6) maser emission toward the central cores in 
NGC6334I. The main NH3(6,6) maser peak is associated with 
the mm peak 2 and elongated along the outflow axis. Therefore, 
we suggest that the strongest line and continuum source in this 
region, mml, is not the driving source of the prominent molec- 
ular outflow, but rather the weaker sub-source, nam2. 



While past observations revealed strong line emission on 
large-spatial scales toward NGC6334I(N), no previous obser- 
vations detected compact gas cores at the center of the region. 
Observations of the highly excited NH3(4,4) to (6,6) lines for 
the first time detect compact warm gas emission from the cen- 
tral main mm peak. The observed line with the highest energy 
above ground, NH3(6,6), is not single-peaked toward that main 
mm continuum core, but it shows a double-horn profile. With 
the current data we cannot differentiate whether this double- 
horn line profile is caused by multiple gas components along 
the line of sight (maybe even a binary system) or whether 
it may be the signature of a potentially underlying massive 
accretion disk. Higher spatial resolution observations of the 
NH3(6,6) or even higher excited lines that can resolve the sub- 
structure of this core are required to solve this question. 
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Fig. 4. NH3(3,3) to (6,6) spectra toward the molecular peak positions associated with the mm continuum peaks mml and mm2. 
The spectral resolution is smoothed to 1 kms"^ and the spatial resolution to 2.8" x 2.2" (the same as the spatial resolution of the 
lower energy lines presented in iBeuther et alj|2005h . The additional strong peaks at positive velocities in the NH3(6,6) spectra 
are a maser discussed in ^3.1.21 The solid lines show fits to the spectra. The (3,3) and (4,4) lines are optically thick so that 
reasonable fits are not possible. The data to this Figure are also available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr ( 130.79.128.5) or via \http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A\ 
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Fig. 5. NH3(3,3) and (6,6) maser emission in NGC6334I is shown in the left and right panel, respectively. Velocity ranges of 
integrated emission are given in each panel and the markers are the same as in Fig. [3] The contour levels are from 15 to 95% 
(step 10%) of the peak emission, which is 112.5 and 126.8 mJ y beam' ^ for the (3,3) and (6,6) maser lines. The dashed contours 
show the Ultracompact Hii region from iBeuther et al.l (l2005h . contoured from 10 to 90% (step 10%) of the peak emission of 
1 197 mJy beam The right panel again shows a smaller region as shown by the inlay in the left panel. The synthesized beams 
are shown at the bottom-right of each panel. The data to this Figure are also available in electronic form at the CDS via 



anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat? J/A+A/^ 
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Fig. 7. Integrated main hyperfine component NH3 inversion line images of the the rmal emission t oward NGC6334I(N) are shown 
in grey-scale with solid contours. The (1,1) and (2,2) images are adapted from Beuther et alj (12005 ). The top-row images are 
shown on larger spatial scales than the bottom-row images, the inlay size for the bottom-row images is shown in the (3,3) panel. 
The contour levels are 58(58)348 mJy/beam and 33(33)198 mJy/beam for the NH3(1,1) and (2,2) lines, respectively. The (3,3) 
to (6,6) integrated images are contoured from 25 to 95% (step 10%) of the peak emission with peak values of 12.9, 11.6, 10.0 
and 10.0 mJy beam \ respectively. The markers are labeled in the top row, SMA mm continuum emissio n is from Hunter e t 
al. (subm.), CH3OH class 11 maser emission from Walsh et al. ( 1 99K ) and th e cm continuum emission from lCarral et al.l (l2002h . 
The mm continuum sources mml, mm2 and mm4 from lHunter et al. I (l2006h are labeled in the bottom-left panel. The large cross 
indicates the two main axis of the two identified molecular outflows. The synthesized beams are shown at the bottom-right of 
each panel. The data to this Figure are also available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr 
( 130. 79.128.5) or via \http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A^ 



